Nonlinear Convection in Reaction-diffusion Equations 
under dynamical boundary conditions 



Gaelle Pincet Mailly'^, Jean-Frangois Raulf^'* 

"LMPA Joseph Liouville FR 2956 CNRS, Universite Lille Nord de France 
50 rue F. Bmsson, B.P. 699, F- 62228 Calais Cedex, France 
(+33)3 21 46 56 50 



Abstract 

We investigate blow-up phenomena for positive solutions of nonlinear reaction- 
diffusion equations including a nonlinear convection term dtu = Au — g{u) ■ 
Vu+f{u) in a bounded domain of under the dissipative dynamical bound- 
ary conditions adtu + d^u = 0. Some conditions on g and / are discussed 
to state if the positive solutions blow up in finite time or not. Moreover, for 
certain classes of nonlinearities, an upper-bound for the blow-up time can be 
derived and the blow-up rate can be determinated. 
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Introduction 

We consider the following nonlinear parabolic problem 

dtU = Au — g{u) ■ Vm + f{u) in f2 for t > 0, 
adtU + d,yU = on do. for t > 0, (1) 

0) = Mo > in Q, 

where ^ : M ^^ M^, / : M ^ M, is a bounded domain of with C^- 
boundary dfl. We denote by ly : dfl i— )■ the outer unit normal vector 
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field, and by dy the outer normal derivative. 

These equations arise in different areas, especially in population growth, 
chemical reactions and heat conduction. For instance, in the case of a heat 
transfer in a medium fi, the first equation dtU = Au — g{u) ■ Vu + f{u) 
is a heat equation including a nonlinear convection term g{u) ■ Vu and a 
nonlinear source /. On the boundary dQ, if a is positive, the dynamical 
boundary conditions describe the fact that a heat wave with the propagation 
speed ^ is sent into the region into an infinitesimal layer near the boundary 
due to the heat flux across the boundary (see and 11|). 

There are various results in the literature about the theory of blow-up for 
semilinear parabolic equations, in particular for reaction-diffusion equations, 
see e.g. jof, 10|, 12| and jsf. In this work, we discuss a problem involving a 
nonlinear convection term. Whereas a Burgers' equation has been studied in 
jsj in the one-dimensional case, we now consider a more general convection 
term and we set in a regular domain of M^. After recalling some qualitative 
properties in Section [H we construct a global upper-solution for Problem 
(II]) in Section |2] and we deduce some conditions on / and g guaranteeing 
global existence of the solutions (Theorem H]). In Section |3l we investigate 
two methods to ensure the blow-up of solutions of Problem ([1]). The first one 
is an eigenfunction method applied to the model problem 

dtU = Au — g{u) ■ Vu + u^ inVt for t > 0, 
adtU + = on dVt for t > 0, (2) 

u{-, 0) = Mo in Q, 

with p > 1 (Theorem [TOj) . We also derive some upper bounds for the blow- 
up time. In the second method, we construct a self-similar lower-solution of 
Problem ([T]) which blows up in finite time. We prove the blow-up of solutions 
of the following problem 

dtU = Au- g{u) ■ Vu + e^'" in H for t > 0, 
adtU + dyU = on dVl for t > 0, (3) 

«(■, 0) = Mo in f2, 

with p > 0. Finally, in Section HI we determine the blow-up rate of the 
solutions of Problem ([2]) in the L°°-norm when approaching the blow-up 
time (Theorem [T6|) . 

Throughout, we shall assume the dissipativity condition 

(T > on 9fi X (0,oo). (4) 
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In order to deal with classical solutions, we always assume that the parame- 
ters in the equations of Problem ([T]) are smooth 



a e C\dQ X (0,oo)) 



(5) 



/ e C\R) , / > 0, 



(6) 



and 



The initial data is continuous, non-trivial and non-negative in Q 



(7) 



Mo € C{Q), Uq ^0, Uq> 0. 



(8) 



Let T = T{a, uq) denote the maximal existence time of the unique maxi- 
mal classical solution of Problem ([T]) 



with the coefficient a in the boundary conditions and the initial data Uq. As 
for the well-posedness and the local existence of the solutions of Problem 
(II]), we refer to Q, and Q. From jcj, since the convection term depends 
linearly on the gradius Vm of the solution, the maximal existence time T is 
the blow-up time of the solution with respect to the L°°-norm: 



1. Qualitative properties 

The aim of this section is to compare the solutions for different parameters 
cr and initial data Uq and to summarize some positivity results on the classical 
solutions of Problem ([T]). Let < ai < (T2 be two coefficients satisfying 
condition vq < uq he two initial data fulfilling hypothesis ([8]) and wq 
a function in Co{Q) with < Wq < Vq. Denote by u„-^, u„^, v and w the 
maximal solutions of the following problems 



e c{Q X [0, T)) n c^'\n X (0, t)) 




= 00 
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dtUa2 = 



Am<^2 - 



= 



in Q ioT t > 0, 
on d^l for t > 0, 
in Q, 



dt.v = Av — g{v) ■ Vv + f{v) in Vt ioi t > 0, 
(T2(9tf + dyV = on dVl for t > 0, 

f (■, 0) = fo in O, 



and 



Sfiu = Aw — (7(1^) ■ Vw + f{w) in for t > 0, 
to = on (9^2 for t > 0, 

w{-, 0) = t^o in VL. 

Let T((Ti,Mo), T((T2,Mo), r(o'2,'?^o) and T{wq) be their respective maximal 
existence times. For the reader convenience, we recall some results stemming 
from the comparison principle Q . 

Theorem 1 ([4*1). Under the aforementioned hypotheses, we have 

T{a2,Uo) < T{a2,Vo) < T{wo) 

and 

< w < f < in X [0, T{a2, uq)) . 
In addition, if Uo G with 

Auq - g{uo) ■ Vmo + /(mo) >0 inVl, (9) 

we have 

T{ai,uo) < r((j2,Mo) 

and 

^ ""o-i in Vlx. [0, T(cri, Mo)) • 



Now, using the maximum principle from [2| , we extend some results obtained 
in in the case of reaction-diffusion to our problem with convection. 

Theorem 2. Assume hypotheses (Hj) - (jS]). Suppose that f > in (0, 00) and 
that a does not depend on time 

(7eC\dQ). (10) 
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Then the solution u of Problem ([T]) satisfies 

■u > in n X (0,r(a-, Mo)), 

dtU>Q m n X [0,T(a,Mo)), 

dtU>Q m Hx (0,T((T,Mo))- 
Moreover, for all C, G (0, T(cr, mq)), there exists d> such that 

dtU > d in Q X [C,,T{a,Uo)). 

Proof. Let r G (0,T(cr, mq)). Since u is C^'^(i7 x [0,r]) and because / and 
g are smooth ((|6]) and ([7])), we can define these constants 

C = sup (7('u) and M = sup (?'(«) ■ Vm — /'(m). 

nx[o,T] nx[o,T] 

First, the comparison principle from j2| applied to Problem ([T]) implies u >Q 
in Q X [0, r] since / > by condition (jG]). Thus we obtain 

dtU > Au - g{u) ■ Vu > Au - C\Vu\ in Vl ioi t > 0, 
adfU + d^u = on dVL for t > 0, 

0) = Mo in ^• 

The strong maximum principle from jif implies 

m := _min u = rninMo , 

nx[o,T] n 

and if this minimum m is attained in f2 x [0,t], u = m in Q x [0, r]. Since 
/ > in (0, oo), the first equation in Problem ([1]) leads to m = 0, and we 
obtain mo = 0, a contradiction with equation ([8]). Hence m > m > in 
Ux (0,r]. 

Then, since the coefficients in the equations of Problem ([1]) are sufficiently 
smooth, classical regularity results in |13j] imply that u G C^'^(fi x [0,r]). 
Thus y = dfU G C^'^(n x [0, r]) and satisfies 

f dty = Ay- g{u) ■ Wy - {g'{u) ■ Wu)y + f'{u)y in for t > 0, 
adty + duy = on dVt for t > 0. 

By continuity, condition ([9]) implies ?/(-,0) > in Vt. The comparison princi- 
ple from [2I implies y > in x [0, r]. In order to apply properly the strong 
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maximum principle, we have to introduce w = ye > 0. By definition of C 
and M, we obtain 



dtw > Aw — g{u) ■ Vw > Aw — C|Vw7| in f2 for t > 0, 
adtw + djyW > on dn for t > 0. 



Again, the strong maximum principle from [2!] implies 



m := _min w = Tcn.nw{-,0) 

qx[o,t] n 



and if this minimum fh is attained in x (0,r], w = m in n x [0,r]. In 
particular, if m = 0, we have dtu = in f2 x [0,r], thus u{-,t) = Uq for all 
t G [0, r]. Hence u attains its minimum in f2 x (0, r], which is impossible due 
to the first part of the proof. Thus w and dtU are positive in i7 x (0, r]. 
Finally, let ^ G (0,r). Because y is continuous and thanks to the previous 
point, there exists d > such that y{-,C,) > d in Q. As y satisfies 



Hence y > d in fl x [^,t]. Note that d depends only on ^, not on r. Without 
this step, we only have y > me~^'^ which may vanish as r — )• T{a, Uq). □ 

An important fact comes from the last statement of the previous theorem. 
For any positive solution u of Problem ([T]), the maximum principle implies 
that for any s G (0, T(cr, uq)), there exists c > such that u{-,s) > c in 
fl. Then, consider the solution m of ([1]) with the constant initial data c and 
the same a in the boundary conditions. Theorem [1] implies u < u. Since c 
satisfies equation (|9]), Theorem [2] leads to dtU > > 0. Thus, u can be big 
enough after a long time (maybe it blows up). So does u, even if uq does not 
satisfy condition Q. 




dty = Ay- g{u) ■ Vy - [g\u) ■ Vu + f\u) ]y 



inVLx r], 
on dVt X [,^, r]. 



the weak maximum principle from [2] implies 



_min y = niiny(-,0) . 

nx[o,T] n 
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2. Global existence 

In this section, we give some conditions on the function g in the convection 
term, which ensure global existence of the solutions of Problem ([1]) for various 
reaction terms /. We use the comparison method from [2|. Thus, we just 
need to find an appropriate upper-solution of Problem ([1]) which does not 
blow up. This is our first lemma. 

Lemma 3. Let a > and K > be two real numbers and let t] G C^([0, oo)) 
with rj' > a"^. For any integer 1 < j < iV, the function U defined inQx [0, oo) 
by 

satisfies 



U{x, t) = K exp (^axj + ri{t)j , 



if 

and if 



dtU > AU - g{U) -VU + f{U) mVlfort>Q, 
adtU + dM>0 on dn for t > 0, 

[/(■,0)>0 inU, 

a9j{(^) > for allLO>0 (11) 

CO 

01 

(y{x, t) > —— for all t > 0. (12) 

Proof. A simple computation of the derivatives of U leads us to 

dtU -AU + g{U) ■ VU = (^r]' - a^^U + agj{U)U in for t > 0. 

Since we assume rj' > a"^, hypothesis ([TT]) implies 

dtU -AU + g{U) ■ VU - f{U) > in ^] x (0, oo). 
Furthermore, on the boundary dVt for t > 0, we have 

adtU + dM = (^ar]'{t) +aiyj{x)^U (13) 



by hypothesis (fT2!) since v is normalized. Finally, f/(x, 0) = i^Texp yaxj + 
r]{0) ) > in i7 is clear. □ 



Remark 1. In the case of the Dirichlet boundary conditions, we can use this 
upper-solution with the special choice r] = (see jH]). But for the dynam- 
ical boundary conditions, we must use a positive t] for the time dependence 
because our solutions are not bounded, see Theorem [2l 

Now we can state the following theorems for a nonlinear reaction term / 
growing as a power of u (Problem ([2])), or as an exponential function (Prob- 
lem ©). 

Theorem 4. Let a be a coefficient fulfilling conditions (jlj), ([5]) and such that 
there exists 6 > with 

inf (J > (5 sup 0" /or t > and (snpa{x,-)) G LJ^^(R^). 

Assume uq satisfies condition (jH]). // there exists an integer I < j < N such 
that 

liniinf^^>0, (14) 
then the solution of Problem is a global solution. 

Proof. In view of Theorem O we can suppose that uq is sufficiently big 
such that there exists C > with 

gj{u) > Cvf~^ in for t > 0. 

For riit) = CS~^ / ( sup a{x, s)) ds + CH, we have rj' > and Equation 

Jo ^ x£dn ^ 
( |T2l) is satisfied. Let ii' be a positive number such that 

K > Mo(a;)e-^^^-''(°) for all x eQ. 

Then by hypotheses (I5l), ([H]) and ( |T2l) . the function U defined in Lemma 
E] is an upper-solution of Problem ([2]) since U{-,0) > uq in Q. Using the 
comparison principle from ^ , the unique solution u of Problem ([1]) satisfies 

< u{x, t) < U{x, i) for all x e and t > , 

thus u can not blow up. □ 
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This theorem holds in particular for a nonlinearity g in the form g{u) = 
{aiu''^, . . . , . . . , aj\fu'^^) with at least one integer j such that aj > 

and Qj > p — 1. A similar result can be derived for Problem ([3]): 

Theorem 5. Under the aforementioned assumptions, the solution of Problem 
is a global solution if the convection term g{u) ■ Vu has (at least) one 
component gj satisfying gj{u) = ajc'^-''^ with aj > and qj > p. 

Proof. Thanks to qj > p, condition f lTT]) is fulfilled because ajc'^^'^ > 
ajC^/u for u sufficiently big. □ 

Remark 2. Condition (fT4l) is optimal for Problem ([2]), see Theorems H] and 
[TOl But it can be improved in some special cases, for example, if the reaction 
term is f{u) = ulnu. Lemma [3] implies that all solutions of Problem ([T]) 
are global if one component gj of g satisfies gj{u) > aj Inu. In fact, in that 
case, every positive solution of ([1]) is global, without any assumption on the 
convection term g. We use this more general argument: 

Theorem 6. Assume that a and Uq satisfy conditions (j4]), ([5]) and ([8]). /// 

is positive in (0, oo) and fulfills 



then the solution of Problem ([T]) is global. 

Proof. In view of the comparison principle, we suppose without loss of 
generality that c <|| Uq ||oo- Consider z the maximal solution of the Cauchy 
problem 



Condition (|T5|) implies that z is a global solution in [0,oo). z is clearly an 
upper-solution of Problem ([1]) and the comparison principle from ^ yields 




(15) 




Its maximal existence time Tj satisfies 




< u{-,t) < z(t) in i7 for t > 0. Thus m is a global solution. 



□ 
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Condition ( IT^ on a allows us to consider fast decaying functions a, but, to 
ensure global existence, it is essential that a does not vanish on the whole 
Indeed let us prove the following blow-up result related to the Neumann 
boundary conditions, for a = on 

Theorem 7. Assume that a = 0, Uq fulfills hypothesis (|H]) and f is positive 
in (0, oo) such that 

r-co ^ 

/ dy < OO for some c > 0. (16) 

Then every positive solution of Problem ([1]) blows up in finite time. 
Proof. Let m be a non-trivial positive solution of 

dtu = Au — g{u) ■ Vu + f{u) in for t > 0, 
dyU = ondVL for t > 0, (17) 

u{-, 0) = Mo in f2. 



Using the maximum principle from |2|, we have u{-,^) > in for ^ > 0. 
Hence, without loss of generality, we suppose Mq > c in fl. Now, consider 
the maximal solution z of the ODE z = f{z) with the initial data z{0) = 
inf{Mo(x) / X G fi}. Condition (fT6l) implies that its maximal existence time 
Tz is finite: 

Jz{o) j\y) 

Since Vz = 0, ^ is a lower solution of Problem f|T7|) . Using the comparison 
principle from j2|, we obtain z{t) < u{-,t) in for t > 0. Thus, u must blow 
up in finite time with < T < T^. □ 

Remark 3. This section illustrates the damping effect of the dissipative 
dynamical boundary conditions: we have shown that for nontrivial a > 
the blow-up time of the solutions of Problem ([T]) can be strictly greater than 
the ones under the Neumann boundary conditions. 



3. Blow-up 

In this section, we investigate the blow-up in finite time for the solutions 
of the Problems (|2]) and ([3]). Let G be a primitive of g and suppose that 
there exist a > and q < p such that 

G{u) < au" for a; > 0. (18) 
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By applying the eigenfunction method (see j^, [9|] and [12|]), we obtain some 
conditions on the initial data Uq which guarantee the finite time blow-up 
and we derive some upper bounds for the blow-up times. Henceforth, we 
denote by A the first eigenvalue of —A in Hq{Q) and by if an eigenfunction 
associated to A satisfying 

^ e H^{n), < V? < 1 in fi. (19) 

Theorem 8. Let a > 0, 1 < q < p , m = p/ {p — q) and suppose G satisfies 
condition ( [T8|) . Assume hypotheses @ - are fulfilled. If 

[ uoip"" dx > {2\n\P-^C)p (20) 
Jq 



C={p~ im [ ^) + ( '-'a^ I iV^r dx , 



with 



then the maximal classical solution u of Problem ([2]) blows up in finite time 
T satisfying 

T < ^-f^ =: T. (21) 

(p-l)(|^]|l-p(^Mo</^-rfxj -2Cj 

Proof. Define 

M{t) = [ wyp™ dx. 



Thus, 

M(t)= [ Auifi'^dx- I g{u) ■Vuifi"'dx+ I u"^"^ dx. 
Jn Jn Jn 

First, we prove that 

J Auip'^dx > -mA|fi|'^(^ J uPifi^'dxy. (22) 

Note that in the case where m G (1,2), yj™ ^ C^(fi). Thus, we consider 
J^{(p + e)"^Au dx, with e > 0. We have 
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Then Green's formula used twice yields 

[ Au{ip+er dx > -m\ [ u^i^+e)"^-^ dx- f (ud^{{^+eD-d^u{^+ef 
Jn Jn J an ^ 

Since f lT9|) implies 

[ ud^i^"^) ds<0 and [ d^u^"" ds = , 

Jdn JdQ 

letting £ — )■ leads to 

/ Auif"" dx > -mX / uif"" dx. (23) 
Jn Jn 

Since ^ < I, J^u(p^dx < j^ui^p dx and by Holder's inequality, fl22|) holds. 
Now, we show that 

~ j g{u)-Vuip"'dx>-ma(^J | V(^r dx) ^ ^^(^^rfa;)'. (24) 

By Green's formula and by definition of G and if, we have 

- / g{u) ■Vuip"'dx = - [ div{Giu))ip"' dx = m f {G{u) ■ V ^)^'^-^ dx . 
Jn Jn J 

Equation (fTSjl and Holder's inequality lead to 

/ {G{u)-V^)^'^-^dx < a I MV™"^|V(/?|c/a; 
Jn Jn 



< aH |V(^|'"da;y'"(^ j u^ip^^ dx " 



Then by the definition of m, ( l24l) is satisfied. 
Henceforth, introduce 



r \ — 

Ci = m\\Q\'^ and G2 = ma(^J |V<^rcixj". 
Then we obtain 

M{t) > j uP ip"" dx ~ Gi(^ f uPip"" dx) " - C2 ( / uPif"" dx^ ' . (25) 
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Set 

1 £ 

pp pp 

El = — — and £2 = q — • 

Ed qJ" 1 

Recall Young's inequality: for a > and e > 0, a = — < 1 for 

e r se^ 

r,s > 1 with + = 1. It yields 



Ci(^ I uPif^'dxy <- I yPyp"^ dx + 



p 

per 



■■=C3 

In the same way, we have 



with 



Then 

M{t) >- [ vP^'^dx-C 
with C = C3 + C4 > 0. By ( |T9i) and Holder's inequality, we obtain that 

Since M is increasing with respect to t, owing to (1201) we have 



and we can conclude that u can not exist globally. To derive an upper bound 
for the blow-up time, we integrate the previous differential inequality between 
and t > 0. We obtain 



-1 

p-i 



M{t) > (^M{oy-p - (p - 1) (^i^r^" - cM{o)-py 

Hence M blows up before T = M{Oy~P{p - 1)-^ (^1\Q\^-p - CM{0)-p^ \ so 
does u. Thus, T <T. □ 
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With a minor change in the proof of Theorem |S1 we can give another criterion 
on the initial data for the blow-up in finite time of the solution of Problem 



Corollary 9. Under the aforementioned hypotheses, if 



j Mo<^™ dx > max|(4mA)^^|fi|,(4C'|l]|P-i)^;| 



(26) 



\v — a/ /o 



with 

2q ^ 
^p-q^ Jn 

then the classical maximal solution of Problem ([2]) blows up in finite time T 
satisfying 

(p - f^UQif"' dxY - 2mX(^ J^Uo(p"'dx^ - 2C"j 

Proof. As in the proof of Theorem [8|, we consider M(t) = uip^ dx satis- 
fying Equation ([25]). We use ([23]) instead of ([22]) 

M > -mXM + j uP^"" cix - C2 ( y u""^"" dx^ ' 

- p 
where Co = ma\ L |V(p|™'(ix ) Proceeding as before with e'^ = 5 — , 

we obtain the differential inequality 

2' ' 

where C is defined as in the assertion. Then 

M > m(^-^^^Mp-^ -m\) + -\n\^-PMP-C' 
- V 4 / 4' ' 

and the blow-up occurs in finite time by ([26]) . To obtain the upper bound for 
the blow-up time, it suffices to integrate the ordinary differential inequality 



M 



> MP - mXM{OY-P + ~ C'M{Oy 



□ 
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We can note that Conditions f l2Up and 0261) on the size of the initial data 
are necessary only to derive an upper bound for the maximal existence time. 
Thanks to Theorem [2], we obtain: 

Theorem 10. Let q < p and suppose G satisfies 

hmsup < oo. 

Assume that a and uq satisfy conditions (jlj), ([5]) and ([8]). All the positive 
solutions of Problem ([2]) blow up in finite time. 

Proof. Let m be a positive solution of Problem ([2]). Theorem [2] permits us 
to ensure that there exist to > and C > such that M(-,to) is big enough 
to satisfy Equation f l20l) and G{u) < Cu'^ in Q ioi t > to. Thus applying 
Theorem [8]to v{x, t) = u{x, t + to), we prove that v blows up in a finite time 
satisfying (1211) . Hence, u blows up in a finite time Tu = t^ + T^. □ 

Now, we present another method to prove the blow-up of positive solutions 
of Problems ([2]) and ([2]) • As in [Tg] , we first construct a lower bound for the 
solutions of the following problem 

dtU = Au — /i| VmI"^ + Ku'P in f2 for t > 0, 
adtU + d^u = on (9fi for t > 0, (27) 

0) = uo in Q, 

with p > q, fi > and k > 0, which blows up in finite time. 

Lemma 11. Let to <ind e be two real numbers with < to < V^- Define in 
VL X [to, l/e) the function V by 

Vix,t) = ^--rWi ^J^^Y 



etf V(l-et)^ 

1 p-q 
-2' q(p-l) 



where k = l/{p - 1), < m < min{|, ^ }, Wiy) = 1 + | - |i /or 



y > with A > and e < yjfxji)' // 1/^ ~ is small enough, and if 
Uo > ^('jto) in Q, then every positive solution of Problem (l27j) blows up in 
finite time. 



Proof. We refer to the proof of Theorem 36.2 from [1^ in the 

1. □ 
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Theorem 12. Under hypotheses (jlD - ([H]); assume that there exists q > 1 
with q < p such that 

limsup^^ < cx). (28) 

Then all the positive solutions of Problem ([2]) blow up in finite time. 

Proof. Let m be a positive solution of (E]). Thanks to Theorem [21 we can 
suppose that uq and u are sufficiently big to verify uq > to) in ^) where V 
is the function defined in Lemma [TTl According to condition (!28|) . there exists 
/3 > such that \g{u)\ < (3u^~^ in VL ioi t > 0. Then, Young's inequality 
implies 

\g{u) ■ Vu\ < /Jm"-^! Vm| < Ku^ + /ilV^r, 

with K = u = —B^ and m = — ^--r. By definition, we have m < p. 
Thus, we obtain 

(9tu >Au- /ilV^r + (1 - k)m*' in for t > 0, 



and ti is an upper-solution of Problem fl271) . and m blows up in finite time 
according to Lemma [TTl D 

Theorem 13. Assume a and Uq satisfy conditions @ - i^. If 

limsup < oo, 

then all the positive solutions of Problem ([3]) blow up in finite time. 

Proof. Let m be a positive solution of ((31) and define v = e'^'^ with 7 G {q,p) 
and 7 > 1/2. As in the previous proof, we suppose that u is sufficiently big 
such that for some C > 

|^(m)| < Ce-^" in for t > 0. (29) 

Computing the derivatives of v, we obtain 

dtv = Av |Vf P — g(u) ■ Vv + 'jv~ in fl for t > 0. 

V 

Using condition ( l29l) . we obtain 

1 I 1 9 £ I P+7 

(9tt' > Av |Vf I — Cvi |Vw| + 7f ^ in for t > 0. 



V 
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Young's inequality 

^ , ^ ,2 1 2. 

Cvy\\/v\< — Vf H — fT, 

leads to 

2 + ,2 £±2 1 22 

dtV > Av 1 Vf I + 7f T V 'I in Q for t > 0, 

since f > 1. Morevover, we have 

p+7 1 22 1 P+7 

7^ 7 y ^ > (7 )V T 

by definition of 7. Thus, in Q for t > we obtain dfV > Av — /i|Vf p + (7 — 
^)v~ with yU = (2 + C^)/2. On the boundary we have 

adtV + dyV = 7f {adtU + j = 0. 

Without loss of generality (see Theorem [2]) , we can suppose that f(-,0) > 
V{-, to) in Q where V is the function in Lemma fTTl Hence, by the comparison 
principle from j2|, v is an upper-solution of Problem f l7r|) and must blow up 
in finite time according to Lemma [11] Thus u blows up in finite time. □ 

Remark 4. In this section, we point out the accelerating effect of the dy- 
namical boundary conditions, in comparison with the Dirichlet boundary 
conditions. Indeed, we prove that, even if the initial data Uq is small, the 
solutions of Problem ([2]) blow up in finite time. But, if we replace the dynam- 
ical boundary conditions by the Dirichlet boundary conditions in the second 
equation of Problem it is well known that the solutions are global and 



decayto if the initial data are small enough, see for instance references [17 
and [iSj and the following theorem: 

Theorem 14. Let g be any non-negative function in C(R, R^) and uq be an 
initial data satisfying ([8]) and uq = on dfl. Then for any p > 1 + 2/N, the 
problem 

g{u) ■ Vm + in Q for t > 0, 

on dn for t > 0, (30) 
in fl, 

admits global solutions if Uq is small enough. 
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Proof. We use this well known upper-solution (see [15| and references 
therein) defined in x [0, oo) by 



U{x,t) = Ait + iyexp 



X 



4(t + l) 

where 7 = and A = ^ — 7 j - In ^5 for t > 0, U satisfies 

d,U-AU + giu).WU-U^=(^^-{:^^-U^^-^]u. 

^ ^ V2(t+1) tr 2(^ + 1) / 

Up to a translation, we can suppose that Xj < in f2 for all 1 < j < N. 
Thus 



dtU -AU + g{u) -VU-U'^y 



\2{t + l) j 



Then by definition of 7, N — 2'j > 0, and ^ < Hence, choosing A 

small enough, we are led to 

dtU -AU + g{u) ■ Vf/ - > in ^] X (0, 00). 

On the boundary, U > is clear. Finally, we just have to consider an initial 
data Mo with uq < U{-,0) in Q, and the comparison principle implies that 
the solution u of (!30!) is bounded from above by U. Thus, u is global and 
goes to as t — )■ 00. □ 

4. Growth Order 

In this section, we are interested in the blow-up rate for Problem ([2]) when 
approaching the blow-up time T. For the convection term, we assume that 

g{u) = {gi{u),---,gn{u)) with gi{u) = Vi = 1, ■ ■ ■ , n, 1 < g G R. (31) 

First, we derive a lower blow-up estimate for p > g -|- 1, valid for any non- 
negative initial data uq € C{Q). 

Lemma 15. Let p > q + 1, and assume hypotheses (jl]) - ([8]). Then the 
classical maximal solution u of Problem ([2]) satisfies 

lk(-,t)||oo> (p-l)^(T-t)^^ 

/or < t < T. 
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Proof. Let t e [0,r). Denote by C ^ ^H(O)'^i)) ^he maximal solution of 
the IVP 

C = in(0,ti) 

C(0) = h(-,t)||oo 

withti = Introduces G C(nx[0,T-t))nC^'^(nx{0,T~t)) 

defined by v{x, s) = u{x, s + t) for x E Q and s G [0, T — t). Then v is the 
maximal solution of the problem 

dtv = Av - g{v) -Vv + vP in for < s < T - t, 
odtv + dyV = on 9i7 for < s < T - t, 

t;(-,0) = inH. 



The comparison principle from [2[ implies that ti < T — t. □ 

This result remains valid for Problem ([T]) as soon as blow-up occurs. We just 
need a positive function / such that an explicit primitive of j is known. 
We follow the technique developed in Theorem 2.3 in Jsj for an one-dimensional 
Burgers' problem and inspired by Friedman & McLeod to prove that 
the growth order of the solution of Problem ([2]) amounts to — — 1) for 
p > 2q + 1 with 4p{p — l){p — 1 — 2q) > Nq"^, when the time t approaches 
the blow-up time T. 

Theorem 16. Suppose conditions (jl]), ([H]), (flOl) and (I3T1) are fulfilled. For 

p>2q + l and p{p-l){p-2q-l) > (32) 

there exists a positive constant C such that the classical maximal solution u 
of Problem ([2]) satisfies 

||<,t)||oo< forte[0,T). 

Proof. Let a > 1 such that 

p{p-\){p-2q-\) = ^^, (33) 
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and choose M > 1 such that 

-2{2q + iy ' 
First, for ^ G (0,T), we shall prove that there exists 5 > such that 

dtu > 5e-^^\uP + 

in f2 X [^,r). Introduce 

J = dtU — 5d{t)k{u) 

with djt) = e~*^* and k{u) = ■u^ + m^'^'*"^. Note that classical regularity results 
from [3 yield J G C^'^ (H x [^,T)). We recall that Theorem [2] implies that 
there exists c > such that dtU > c > in Q x [C,,T). Thus, we can choose 
6 > sufficiently small such that 

J(-,0 > in n. 
J fulfills the boundary condition 

adtJ+d^J = dt{adtu+d^u)-6dk'{u){adtu+d^u)-a6d'k{u) = a6Me-^^'k{u) > 0. 
Furthermore, J satisfies 

dtJ -AJ + g{u) ■ V J - {pu''-^ - g'{u) ■Vu)J = 6dH{u) in H x [^, T), 
where 

H{u) := pvF-^k{u) - k!{u)vF + k!\u)\Vu\^ - -k{u) - k{u)g\u) ■ Vu. 



d 



To prove that H{u) > 0, we shall show that 



qVNu'^'^\Vu\{uP + u'^''^^) < M(mP + + - 2g - (34) 

+ {p{p - + 2g(2g + l)u'^'^~^)\Vu\'^. 

Inequality is trivial in the case where M > q\/Nu'^^^\'Vu\. Now, sup- 
pose that M < qy/Nu'^-'^\Vu\. When qy/Nu'^^^ < 2g(2g + 1)|Vm|, we 
have q^/Nu'^-^u'P\Vu\ < p{p - 1)u'p-^\Vu\'^ and q^/Nu^'^\Vu\ < 2q{2q + 
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^iVwp since p > 3 then fl34p follows. In the case where qy/Nu'^'^^ > 
2q{2q + l)|V-u|, since 

we obtain 



Nq 



u 



Moreover, (1551) yields 

aViVgu^+^lVMl < 2y/p{p- l){p-2q- l)u'^^^\Vu 

< (^Vj9-2g-lM'^+^ - v/p(p - 1) I V 

+2^j9(p- l)(p-2g- 1)m''+^|Vm| 

< (p-2g- 1)|Vm|2. 

Thus, multiplying by we are led to 

aVNqu'^-^\Vu\uP < (p - 2g - 1)^^+^^ + p(p - l)^^-^! Vm^ 

and by (l35l) . the inequality ( l34l) holds. Finally, we can conclude by the 
comparison principle from [2[ that J > in fix T), in particular, dtU > eu^ 
with £ > 0. 

Now, we shall derive the upper blow-up rate estimate of ||M(-,t)||oo for t G 
[^,T). For each x G fi, the integral 



uP{x, s) J^^^^t) r]P 

converges as r — )> T. Integrating the inequality dtU > eu^ leads to 

eiT — t) < < . 

^ ^ - 1-p -p-1 

-1 

Letting t ^ T implies u{x, t) < l)(T-t))'"' and we can conclude 

as in the proof of Theorem 2.3 from □ 



The latest result can be_improved thanks to Theorem 38.1 developed by 
Quittner and Souplet in [14] in the Dirichlet case. We extend their theorem to 
the case of the dynamical boundary conditions. For the reader's convenience, 
we precise the minor changes appearing in the proof. 
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Theorem 17. Assume a > on the boundary dQ and conditions ([H]), f llUp 

and 031 p are fulfilled, and suppose that 

f oo zfN = l 

l<p<Ps:=^N^ ./iV>l 

Assume in addition that h : ]R_|_ x R" — )■ R satisfies the growth assumption 

\h{u,o\<Coii + \ur + \^n, (36) 

for some 1 < Pi < p and 1 < p2 < '2p/{p + 1). Then, any nonnegative 
classical solution of 

dtu = Au + h{u, Vu) + uP m n, < t < T, , . 

_ adtu -d^u = on dn,0 <t <T ^ > 

satisfies 

u{x,t) + \Vu{x,t)\^ < C{l + t^ + (T-t)^), (38) 
for xEfl,0<t<T, with C > depending on p,pi,p2, Co, ^■ 

Proof. We assume that the estimate (1380 fails. Then as in Theorem 38.1 
from [14], there exist sequences Uk, G (0, oo), Xk E Q and G (0,Tfc) such 
that Uk is a solution of ( 1371) in (0, T^) with: 

Mk = u— , Mk{xk,tk) > 2k , Mk{xk,tk) > 2k min{sfc, Tk — tk} ^ 

2 

and Afc = 1/Mk{xk,tk). The function Vk{y,s) := \lr^Uk{xk + XkV.tk + A^s) 
defined in IDk = (^\k\^ - ^k) n {\y\ < k/2}^ x {-k^/A, k^/A) satisfies 

dgVk = Avk + hk + vl in Dk, 

f-^dsVk + dyVk = iiye Xk\dn - Xk), \y\ < |, \s\ < f . 
with 

2p -2 -(p+1) 

/ifc(y,s):=A-^/i(Ar'^^fe(2/,s),A,''-^ Vvkiy,s)) , 
^^r^(0) + |Vt;fc|^(0) = l, 

and 

Vk"" + |Vt;fc|f+i < 2 in Dfc . 
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Condition implies \hk\ < CA^ with m = min{^^^5fi^, ^^^^}. Then, 
as in the proof of Souplet and Quittner, we obtain a subsequence of (vk) 
converging in C^+^'^/^j^]^^ x M) with < a < 1 to a nonnegative function 
V. Observe that in the hmit case — )■ 0, the dynamical condition dgVk + 
^dyVk = on the time lateral boundary of Dk, implies that v does not 
depend on time t at space infinity, that is f = f(y,0) > 0. Two cases can 
occur. In the first case f is a solution of 

dtu = Au + vF for a; G M^, t > 0, 

and in the second one, v satisfies 

dsV = AyV + lyP y G, He, s G M, 
v>0 ye dHc, seR, 

with He := {y e / yi > -c}. Since w^(0) + \Vv\^{0) = 1, is 
non-trivial and v and Vv are bounded. Hence, we are led to a contradiction 
with non-existence results (see Theorems 21.1 and 21.8 from [mI]) for positive 
solutions of dfU = Au + in unbounded domains. □ 

Thanks to this last theorem, we can improve the result obtained in Theorem 
[TBI Indeed, we derive a blow-up rate estimate for p > 1 -|- 2g > 3 and p < Pb 
where the constant C involved in the upper bound is independent on the 
initial data unlike in Theorem [T6l 

Theorem 18. Under the aforementioned hypotheses, suppose p > 1, q > 1 
and 

2q + 1 < p < Pb ■ 
Then, the maximal classical solution of Problem ([2]) fulfills 



2 



sup + |Vm|~ j <Cyl + t~ + (T - t)~ 

for t e [0, T), with a constant C such that < C = C{p,q,fl). 
Proof. Young's inequality yields 



\9{u) ■ ^1 < ^/Nu%\ < [j-y)"^ + ^1^1 



P2 



for C e M^. Equation (IMD is satisfied with pi = 2±f±i and p2 = 

thus Theorem [T71 permits to conclude. □ 
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Since p > 1 + 2g > 3, this theorem apphes when ps > 3, i.e. for = 1, 2 
and 3. In these cases, our last result completes Theorem (TH] giving a blow-up 
estimate for p G (3,pb). In particular for q = I, Equation fl32|) is fulfilled for 
any p > ps, then the blow-up rate for the solutions of Problem ([2]) is fully 
determinated for all p > 3. 
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